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ABSTRACT 
The  elastic  transient  response  of  an  impulsively  loaded  circular 
ring  is  examined.   The  numerical  stability  criterion  associated  with 
a  finite  difference  formulation  of  the  problem  is  briefly  discussed, 
and  a  method  to  reduce  the  stability  restriction  is  presented.   This 
method  assumes  that  the  density  of  the  ring  is  greater  in  the  tangen- 
tial direction  than  it  is  in  the  radial  direction.   The  analytical 
solution  to  the  ring  equations  with  the  artificial  tangential  inertia 
is  determined  for  a  specific  impulsive  load.   This  solution  is 
compared  with  the  solution  to  the  ring  equations  with  the  correct 
tangential  inertia.   The  comparison  shows  that  on  the  basis  of  pre- 
dicting peak  values  of  stress,  the  results  for  the  ring  with  the 
artificial  tangential  inertia  are  very  close  to  those  of  the  unmod- 
ified ring.   Thus,  the  method  reduces  the  numerical  stability  require- 
ments and  does  not  significantly  affect  the  principal  results. 
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CHAPTER  I 
INTRODUCTION 

The  finite  difference  method  of  analysis  has  long  been 
recognized  as  an  effective  and  practical  means  for  solving  partial 
differential  equations  on  a  digital  computer.   The  accuracy  asso- 
ciated with  this  method  can  be  adjusted  by  changing  the  incremental 
sizes  of  the  independent  variables  involved.   In  the  limit,  when 
the  sizes  of  the  increments  approach  zero,  the  solution  to  the 
finite  difference  equations  becomes  identical  to  the  solution  to 
the  corresponding  differential  equations.   Consequently,  the  engineer 
may  improve  the  accuracy  of  his  finite  difference  solution  by  choos- 
ing sufficiently  small  values  for  the  incremental  sizes  of  the  inde- 
pendent variables,  with  greater  accuracy  accompanied  by  increased 
execution  time  on  the  digital  computer.   However,  there  are  many 
problems  where  the  finite  difference  method  imposes  a  restriction 
on  the  sizes  of  the  increments.   This  restriction  is  known  as  the 
stability  criterion  and  arises  from  a  purely  numerical  phenomenon. 
The  stability  criterion  defines  a  critical  limit  for  the  ratio  of 
the  increment  sizes  above  which  a  numerically  stable,  valid  solu- 
tion will  result.   A  value  of  the  ratio  below  this  limit  will  lead 
to  a  divergent,  meaningless  solution. 

The  stability  criterion  for  the  finite  difference  dynamic 
analysis  of  the  thin,  homogeneous,  elastic,  circular  ring  has  been 
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defined  by  Ball.    It  can  be  simplified  to 


where  A  is  the  radius  of  the  ring;  h  is  the  thickness  of  the  ring; 
D  =  Eh,  the  stretch  rigidity;  ^      is  the  ring  density;  and,  A  (t) 
and    ^7"   ^i^s  the  angular  and  time  increments,  respectively.   A 
meaningful  numerical  solution  for  the  transient  response  of  the  rin^ 
is  possible  only  when  the  above  restriction  on   Z^(p        and   AT 
is  observed.   A  physical  interpretation  of  the  stability  criterion 
says  that  the  finite  difference  solution  must  propagate  around  the 
ring  at  least  as  rapidly  as  a  stress  wave. 

Once  the  size  of   A^   is  chosen,  the  maximum  size  for  the 
time  increment  becomes 


(ATU    =    (A^)A^^ 


An  accurate  solution  for  the  ring  response  requires  a  small  value 
for  the  angular  increment,   A^  •   Therefore,  according  to  Equation 
(1)»  (  AT   )     is  automatically  restricted  to  a  small  value;  and 
hence,  large  computation  times  will  be  required  on  the  computer. 

To  reduce  the  restrictions  on  the  numerical  stability  of  the 
ring  response  a  method  is  proposed  whereby  the  maximum  time  increment 
may  be  increased  without  producing  an  unstable  solution.   According 
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to  Equation  (1),  (  AT  )^au   is  determined  once  A^*  is  chosen. 
However,  if  the  density  ^   is  considered  to  be  a  variable  in  the 
stability  criterion  then  (  AT  )f/»flx  ^"^^^   ^^   proportional  to  the 
square  root  of  the  density  at  a  particular  value  for  A  ^  •   The 
density  in  Equation  (1)  is  a  measure  of  the  ring's  inertia  in  the 
direction  of  the  propagation  of  stress  disturbances  around  the  ring, 
that  is,  the  tangential  direction.   The  method  proposed  here  modifies 
the  ring  by  introducing  an  artificial  tangential  density  that  is 
larger  than  the  actual  ring  density.   Thus,  the  ring  inertia  will  be 
characterized  by  the  density  ^    in  the  radial  direction  and  by  a 
different  density,  X  p      >  in  the  tangential  direction.   At  a  given 
value  of   A  db  >  the  maximum  time  increment  for  a  stable  response 
of  the  modified  ring  will  be  proportional  to  the  square  root  of  ^ P 
and,  hence,  will  be  larger  than  the  critical  time  increment  for  the 
original  ring  by  the  factor  y  ^  . 

For  example,  for   ^  =  4  the  time  increment  for  stable  response 
of  the  modified  ring  is  twice  as  large  as  the  time  increment  for  the 
original  ring  for  the  same   A®  •   Thus,  the  computation  time  can 
be  cut  in  half.   The  decreased  computational  time  afforded  by  the 
proposed  method  is  a  significant  improvement  on  the  conventional 
finite  difference  method  of  analysis.   However,  the  introduction 
of  the  artificial  tangential  inertia  will  also  lead  to  an  incorrect 
solution  for  the  ring  response.   Thus,  the  value  of  the  method  will 
depend  upon  the  accuracy  of  the  solution  of  the  modified  ring. 
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This  paper  investigates  the  effect  of  an  artificial  tangential 
inertia  on  the  analytical  solution  for  a  particular  ring  problem. 
Solutions  are  obtained  for  the  elastic  response  of  an  impulsively 
loaded  circular  ring  with  and  without  artificial  tangential  inertia. 
The  maximum  stress  in  the  ring  is  computed  for  the  original  ring 
(  Jf  ■  1)  and  for  rings  with   ^  "2,  3,  and  4.   A  comparison  of  the 
maximum  stress  for  each  ring  demonstrates  the  value  of  the  method. 

This  study  was  completed  during  the  1967-1968  academic  year  at 
the  Naval  Postgraduate  School,  Monterey,  California.   Acknowledge- 
ment is  gratefully  made  to  Professor  Robert  E.  Ball  of  the  Aeronautics 
Department  for  his  guidance  and  consultation. 
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CHAPTER  II 
A  DISCUSSION  OF  THE  RING  EQUATIONS 


This  thesis  considers  a  thin,  circular  ring  composed  of  a 

homogeneous  material.   The  governing  equations  are  derived  in 

Reference  4  and  are  repeated  here  for  completeness.   Figures  1  and 

2  illustrate  the  general  geometrical  characteristics.   The  radius  A 

is  measured  to  the  middle  surface  of  the  ring,  which  lies  at  a  dis- 

u 
tance  —-  from  both  edges.   A  differential  element  on  the  ring,  O  (p     , 

is  located  at  an  angle  ^  from  the  bottom  of  the  ring,  where  <^     is 
measured  positive  counter-clockwise.   The  displacement  of  a  point 
on  the  middle  surface  is  described  by  a  radial  and  a  tangential 
component,  W  and  V,  respectively.   The  radial  distance  from  the 
middle  surface  to  some  arbitrary  point  on  the  ring  is  given  by  Z  in 
Figure  2. 

A  two-dimensional  problem  will  be  assumed  with  regard  to  the 
internal  stresses  of  the  ring.   Figure  3  illustrates  the  forces 
acting  on  the  element   Q^  ^^^    to  the  normal  and  shear  stress  dis- 
tributions on  the  cross-sections.   The  forces  N,  M,  and  Q  are, 
respectively  J  the  normal  force,  the  bending  moment,  and  the  radial 
shear  force.   These  forces  are  stress-resultants  referred  to  the 
middle  surface  of  the  element.   The  definitions  of  N  and  M  are  given 
by  the  following  expressions 


,\ 


N  = 


fl 


c^t  dz 


'i> 


2. 


£  €^dz. 


(2a) 


2 
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FIG.  I 
GENERAL  RING  GEOMETRY 
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FIG.    2 
RING  MIDDLE  SURFACE 


Middle 
Surface 


FIG.   3 
TWO-DIMENSIONAL  STRESS-RESULTANTS 


17 


h 
IT 


/I 


M  = 
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Ee^  2dz 
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(2b) 


Z 
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where  <J^      and   £i  are  the  stress  and  strain  on  the  cross  section 
at  a  distance  Z  from  the  middle  surface. 

In  order  to  obtain  an  expression  for  the  strain  £^  in  terms 
of  Z,  two  additional  assumptions  are  made: 

1.  A  line  normal  to  the  middle  surface  before  deformation 
will  be  normal  to  this  surface  after  deformation. 

2.  A  point  on  the  normal  to  the  middle  surface  will  have  a 
constant  distance  Z  from  this  surface,  both  before  and 
after  deformation. 

Accordingly, 


_   /    ^V 


^-^  =  ^ 


4- 


#w 


w 


*      A    a<^       /IfAfZ)    ^c^       AfZ 


(3) 


Substituting  Equation  (3)  into  the  definitions  for  N  and 
M,  Equations  (2a)  and  (2b)  ,  leads  to 


» •  x(ij  - ")  -  f  (1^ ' " 


M  =  - 


j1 


( 


■f  W 


) 


(4) 


(5) 
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where  D  =  £/?  /f  -  ^ 

The  linear  equilibrium  equations  for  the  element  follow  directly 
from  an  examination  of  Figure  3.   They  are 


^   -Q    -    fAh  —  (ea) 

where  rotary  inertia  is  neglected. 

Eliminating  the  radial  shear  Q  from  these  equations  leads  to 

-M  >  _L_  _^M_  -    ^Ah^ 

a^      A    <^4)      ^      ^T^  (7a) 

Substituting  Equations  (4)  and  (5)  ,  relating  the  stress^resultants 
to  the  displacements  W  and  V,  into  Equations  (7a)  and  (7b),  gives 

These  are  the  differential  equations  that  describe  the  free 
vibrations  of  the  thin,  two-dimensional,  elastic  ring. 
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CHAPTER  III 
SOLUTION  FOR  THE  IMPULSIVELY  LOADED  RING 

In  order  to  adequately  test  the  proposed  method  for  increasing 
the  critical  time  step,  an  impulsive  load  that  causes  both  large 
bending  and  membrane  stresses  should  be  employed.   Consequently, 
the  impulsive  load  illustrated  in  Figure  4  is  taken  as  the  loading 
on  the  ring  at  time  t  =  0.   The  impulse  is  applied  to  the  top  half 
of  the  ring  and  is  characterized  by  a  cosine  distribution  directed 
radially  inward. 

The  initial  conditions  for  the  radial  cosine  impulse  are  given 
by 

W(<^>  0)=  0  V(<^>0)   =  0 

=  0  -f«<^.<f     (9) 

where  V,    =  -^ 

Since  the  ring  is  unrestrained,  a  net  momentum  will  be  imparted 
by  the  radial  velocity  acting  on  the  upper  half  of  the  ring.   This 
will  cause  a  rigid  body  translation  of  the  ring,  and  hence,  unbounded 
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FIG.  4 
IMPULSE  LOADING  ON  THE   RING 
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values  for  the  displacements  V  and  W.   To  avoid  this  situation,  the 
initial  conditions  must  be  modified  in  such  a  way  that  the  rigid 
body  motion  is  eliminated,  while  preserving  the  internal  effect  of 
the  impulse.   Reference  2  defines  the  following  modified  set  of 
initial  conditions  for  this  loading 

w{<P',o)-o  v(4>,o)'0 


^(*.0)=  -^SIN4> 


-fc0S<p  -f**«f    (10) 

With  this  set  of  initial  conditions,  the  ring  will  remain  fixed  in 
space;  and  hence,  the  displacements  V  and  W  will  measure  actual 
deformations  of  the  ring. 

The  following  solutions  for  W  and  V  are  assumed 


oO 


W(4>.r)=  )    A^(r)  cos  ncf> 


(lla) 


oO 


V{<PjT)   =>       Bnir)    SIN    r\<j>  (lib) 


(12) 


*a   and  b   are  arbitrary  constants, 
n      n 
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Based  upon  these  solutions  and  the  given  initial  conditions  the 
expressions  for  the  bending  moment  and  normal  force  are* 


2  ^ 1         %+/ 


n=2,4,.. 


f'-n^ 


^. 


COS  n  4> 


SIN^r^Od^T 


(13a) 


^o  "  ^-^^f^  =  'S//V  aA>7  -^  ^-^   SIN  {2  60,7    COS  4> 


+ 


C-/)'^^' 


n=2,4,... 


.    (nM)(^rf-Sn') 


-f^  5//y-/^„  6Jor 


■^  <^y)  S/W  9n  <^or 


COS  K)  <^  -  «>  Mc 


(13b) 


where 


o(  ^ 


K 


DA 


"-ta 


Mq  and  Nq  are  non-dimensional  representations  of  the  bending  moment 
and  the  normal  force,  both  of  which  may  be  solved  without  specifying 
the  impulse  magnitude,  1^.   The  other  terms  will  be  described  later. 
Equations  (13a)  and(13b)  characterize  the  dynamic  response  of  the 
elastic  ring  due  to  the  specified  impulse.   Before  proceeding  further, 
however,  a  physical  interpretation  will  be  givaxto  these  expressions. 


*These  expressions  are  derived  in  Appendix  A.   Since  this  is 
the  original  ring,  ^  =  1. 
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Consider,  first,  the  various  mode  shapes  excited  by  the  impulsive 
loading.   Figure  5  illustrates  the  lower  modes  of  vibration  tor  the 
ring.   The  lowest  mode,  n  «  0,  is  characterized  by  a  radial  displace- 
ment only,  where  W  is  the  same  for  all  points.   A  small  change  In 
curvature  in  addition  to  the  middle  surface  straining  will  result; 
and  hence,  N^  and  M,  will  both  have  contributions  from  this  mode  as 
can  be  seen  from  Equations  (4)  and  (5).   The  next  mode,  n  ■  1, 
involves  middle  surface  straining  only.   The  shape  of  the  ring 
remains  the  same,  with  no  change  in  the  effective  curvature  of  the 
ring.   Hence,  this  mode  will  contribute  only  to  the  normal  force, 
No  •   For  n  >  1,  only  the  even  mode  shapes  appear  in  the  expressions 
for  the  stress-resultants.   This  situation  arises  from  the  modified 
initial  conditions  imposed  on  the  ring.   At  time  t  ^  0,  positive 
radial  velocities  are  required  at  the  top  and  bottom  of  the  ring, 
^  *  7r  and  0,  while  zero  radial  velocities  are  prescribed  at  the 


lides,  4**  ^  9  •   Referring  to  Figure  5,  it  can  be  seen  that  only 


2L 
2 

the  even  modes  can  account  for  these  velocities  and  hence,  are 
the  only  ones  excited  by  the  impulse. 

Finally,  examination  of  Equations  (13a)  and  (13b)  reveals  that 
two  frequencies,  f  cja  ^^-^  8  COo  '  ^*"^  associated  with  each  mode 
of  vibration  of  the  ring.   These  frequencies  arise  from  the  fact  that 
the  ring  element,  cl  <b  >   may  displace  independently  in  two  separate 
directions.   From  Appendix  A,  the  frequency  coefficients,  f  and 
g  ,  are  approximately  given  by 
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n=0 


n  =  I 


n=  2 


n  =  4 


FIG.    5 
LOWER  MODE  SHAPES 
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(14) 


Since  a  thin  ring  is  assumed,  o<>  is  very  small,  and  f  <<  g  .   The 

n       n 

frequency  coefficients,  f  and  g  ,  characterize,  respectively,  the 

radial  and  tangential  vibrations  of  the  ring.   Due  to  the  relative 

magnitudes  of  f   and  g  ,  the  expressions  for  M^  and  N_  will  each  be 
n      n 

dominated  by  one  of  the  two  sine  terms  of  the  series.   Hence,  M- 


will  vary  at  a  frequency  given  approximately  by  f  CO©    while 
will  be  characterized  by  a  much  larger  frequency  g  (^    , 


N, 


To  obtain  numerical  data  from  Equations  (13a)  and  (13b) ,  it  is 

necessary  to  select  physical  data  describing  the  ring.   Figure  6 

outlines  the  physical  characteristics  for  a  thin  ring  composed  of 

aluminum.   The  frequency  coefficients,  f  and  g  ,  have  been  calcu- 

n      n 

lated  on  a  digital  computer  at  20  micro-second  intervals  over  a 
period  of  time  corresponding  to  one  complete  cycle  for  M^  .   Stresses 
on  the  inner  and  outer  edges  of  the  ring  were  also  calculated  from 
the  following  non-dimensional  relationships 


(outer   edge) 


£VL 


N.  .  ^  M» 


(15a) 


(inner  edge) 


0^.    - 


_  AcdoTT 


£V. 


CT/     = 


(15b) 
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h  =  .25  in. 


CUq    =  16980  rod/sec. 
D       =  2.74  X  lO'  Ibf./in. 
E       '  10.9  X  10*  Ibf./in.* 


K    =  1.43   xlO'    Ibf.-in. 


P 

a 


3.14  xld*  ^-^ 
3.54  X  10"* 


FIG.    6 
PHYSICAL  RING  DATA 
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The  maximum  values  for  M^  ,  N^  ,  and  C^  found  over  the  time  sampled 
have  been  listed  in  Table  T.* 

Figures  7  through  10  illustrate  the  dynamic  response  of  the 
ring  through  the  first  half  cycle  of  M©  .**  It  should  be  noted  from 
Table  I  that  the  response  at  <p  ^    -^   i  ""^^    ^^  °^  secondary  importance 
and  has  been  omitted  from  the  curves.   Also,  the  response  at  ^  =  TT 
is  identical  in  character  to  that  at  <^  =  0     and  has  been  omitted. 
The  symmetry  of  the  even  modes  make  it  unnecessary  to  determine 
solutions  for  ring  locations  in  the  range  7T  ^  ^  $  27r. 


*The  maximum  (J^    was  observed  to  be  within  a  few  per  cent  of 
maximum  0"^  found. 

**A  positive  M^  indicates  a  negative  bending  moment. 
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TABLE    I 
MAXIMUM  VALUES   FOR  THE    STRESS-RESULTANTS 


TT 
4 

7r 
2 

4 


0  2.22  -129.6  -3.45 


1.80  26.5  1.99 


1.31  125.8  2.57 


1.77  26.5  1.80 


TT  2.39  -129.6  -3.27 
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CHAPTER  IV 

SOLUTIONS  FOR  THE  IMPULSIVELY  LOADED 

MODIFIED  RING 

The  assumption  that  the  ring's  inertia  characteristics  are 
variable  with  regard  to  tangential  vibrations  will  now  be  considered. 
The  ring  equations  must  be  modified  to  account  for  this  assumption, 
and  modified  solutions  analagous  to  the  original  ring  are  obtained 
in  Appendix  A  for  the  impulse  loading  described  in  Chapter  II- 

Figure  11  illustrates  the  assumed  inertia  forces  associated 
with  the  ring  element,  O  (i>     .      A  simple  adjustment  of  the  ring 
equations  will  account  for  this  situation.   The  following  expres- 
sions for  the  stress-resultants  J  derived  in  Appendix  A,  are  obtained 
for  the  modified  ring 


<aO 


M.  =  SIN  od,r  -h 


Y)^2,4,. 


SlN-fn^T 


-2 


Cos  ni> 


COS  i> 


oo 


+ 


2M) 


^+1 


^^^^.,,(n^-i)(^-g„') 


-^^S/Z^fyjCOoT 


-h  <Bn  Sthl  9r}<^T 


COS  n<P  -  ^  /W« 


COD 


34 


FIG.    II 
ASSUMED  INERTIA  FORCES 
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To  determine  the  effect  of  the  factor  ^  on  M^^  and  N^,  ,  consider 

first  the  change  in  the  frequency  coefficients,  f  and  e"  .   As  deter- 

■^  n     °n 

mined  in  Appendix  A,  these  coefficients  are 


fn    =    n(n^-/)  \j-^  (18) 


9n    -  V^  (^^> 


For  this  analysis,  it  will  be  assumed  that  JT  takes  on  values  between 
1  and  4.   Therefore,  f  will  not  be  significantly  affected  by  6    , 
especially  for  the  higher  mode  shapes.   However,  "g     will  be  very 
sensitive  to  a  change  in  ^  ,  as  can  be  seen  in  Equation  (19).   Figures 
12  and  13  graphically  illustrate  the  error  in  the  modified  frequency 
coefficients  as  a  function  of  ^  .   A  physical  interpretation  of  the 
factor  y  will  verify  the  above  results.   As  0  increases,  a  progres- 
sively heavier  system  is  seen  in  the  tangential  direction,  giving 
rise  to  a  response  of  lower  frequency,  g"  COq  •   The  radial  vibrations, 
characterized  by  f  Ol)©  >    should  not  be  significantly  altered  by  ^ 
since  the  density  in  this  direction  is  always  given  by  the  original 
ring  density,  P  . 

Values  for  the  modified  frequency  coefficients  have  been  calcu- 
lated and  listed  in  Appendix  B  for  the  ring  defined  in  Figure  6. 
Again,  it  can  be  seen  that  f  <  <  g"  .   M^  and  ITo  ,  then,  will  be 
approximately  described  by  sinusoidals  at  frequencies  f  COo  and  g"  cO  > 
respectively. 
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Examination  of  the  expressions  for  m'^,  and.  W^   reveals  the 
dominant  terms.   For  the  bending  moment,  Equation  (J.6),  the  lower 
frequency  sine  term  will  dominate.   The  magnitude  of  the  coeffi- 
cient for  this  term  will  decrease  as  the  product,  ^g'^,    increases, 
corresponding  to  increasing  values  for  ^  .   Thus,  the  time 
response  for  the  bending  moment  will  decrease  in  magnitude  as  )f 
increases.   The  wave  form  for  the  response,  however,  will  not  change 
significantly.   Examining  Equation  (17)  in  a  similar  manner,  it  can 
be  seen  that  both  the  wave  form  and  the  magnitude  of  the  time  response 
for  No  will  change  as  ^  increases. 

Figures  14  through  19  illustrate  the  curves  for  M.^  ,  N^  ,  and 
OL  for  three  values  of  ^at  <p  "  0    •   Figures  20  through  22  sho\ 


)W 


M"o   for     4>    '■~.*     The   values   obtained   for  M^  and  N^  at    4>   «  ■—-  and 

Cm  J^Lf 

^  *  -2r- were  found  again  to  be  small  compared  with  other  locations 
on  the  ring  and  have  been  omitted;  Also,  the  curves  for  cp  *  7T 
have  not  been  included  because  of  the  similarity  with  the  response 
for  45-0, 

Tables  II  through  IV  list  the  maximum  values  observed  for   the 
solutions  o£  M"g,  ,  rJ^  and  (JJj  .   The  corresponding  values  for  the 
unmodified  ring  have  also  been  included  for  comparison. 


*!T<j  and  CT^'  at  <^  *  -^  are  not  included  because  of  their  small 
magnitudes. 
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TABLE  II 
MAXIMUM  BENDING  MOMENTS 


c^    Jf=i    }r=2   2r=3    ^=4 


0  -129.6        -129.5      -124.9  -117.7 

125.8  118.3        110.4  103.2 

71"  -129.6        -129.5      -124.9  -117.7 

TABLE  III 
MAXIMUM  NORMAL  FORCE  (MAGNITUDES) 


^ 

y-/ 

y=2 

3r=3 

y=^ 

0 

2.22 

2.38 

2.50 

2.48 

TX 
2 

1.31 

1.48 

1.14 

1.13 

n 

2.39 

2.43 

2.46 

2.48 

TABLE 

:  IV 

MAXIMUM  OUTER  EDGE 

STRESS 

4> 

jr»; 

}r^2 

I 

=  3 

5=4 

0 

-3.45 

-3.49 

-3. 

28 

-3.57 

2 

2.57 

2.38 

2. 

20 

2.18 

TV 

-3.27 

-3.09 

-3. 

31 

-2.98 
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CHAPTER  V 
CONCLUSIONS 

A  comparison  of  the  response  of  the  original  ring  with   the 
results  for  the  modified  ring  with  ^  =  2,  3,  and  4  can  be  made  using 
the  data  presented  in  Tables  I  through  IV  and  the  response  curves 
shown  in  Figures  14  through  22.   A  study  of  the  time  response  for 
the  bending  moment,  M^  and  Mq  ,  reveals  that  the  solutions  to  the 
modified  ring,  in  general,  range  from  fair  (  ^  ~  2)  to  poor  (  J  =  4) 
approximations  to  the  solution  for  the  original  ring.   Average  errors 
range  from  10  per  cent  for  JT  =  2 ,  to  about  40  per  cent,  for  S*  =  4. 
However,  the  percentage  error  associated  with,  each  modified  curve 
is  nearly  constant,  and  the  wave  form  is  identical  to  the  unmodified 
curve.   It  would  appear,  then,  that  the  data  from  a  modified  solu- 
tion could  be  increased  by  some  constant  factor  to  yield  a  very 
good  approximation  over  the  first  quarter  cycle. 

Based  upon  the  above  comparison,  the  proposed  method  appears  to 
introduce  unacceptable  errors  into  the  solution  for  J"  >  2 .   However, 
the  actual  time  response  for  the  bending  moment  is  of  little  signifi- 
cance in  a  dynamic  analysis,  since  the  maximum  value  of  bending 
moment  is  usually  the  desired  result.   The  maximum  bending  moments 
observed  for  all  solutions  have  been  listed  in  Table  II.   Based  on 
these  values  the  solutions  for  the  modified  ring  are  very  good 
approximations  to  the  results  for  the  unjnodified  ring.   Errors  in  the 
maximum  bending  moment  of  about  4  per  cent  and  8  per  cent  can  be 

seen  for  ^  values  of  3  and  4,  with  essentially  no  error  observed 
at  y  =  2. 
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The  effect  of  a  heavier  system  in  the  tangential  direction 
is  readily  observed  in  an  examination  of  the  various  normal  force 
response  curves.   Both  the  wave  forms  and  the  maximum  values  are 
affected  by  a  change  in  the  factor  0  .   As  anticipated,  the  tangen- 
tial vibrations  decrease  in  frequency  as  J  increases.   This  can  be 
seen  from  the  spacing  of  the  beats.   Despite  the  obvious  changes 
in  the  normal  force  curves  for  the  modified  ring,  the  maximum  peak 
values,  listed  in  Table  III,  are  not  significantly  altered  as  S 
changes . 

The  primary  concern  of  any  dynamic  analysis  is  usually  the 
maximum  stress  that  occurs  in  the  ring.   Since  both  the  bending 
moment  and  the  normal  force  contribute  to  the  stresses,  the 
maximum  value  will  be  found  on  one  edge  of  the  ring.   The  relative 
importance  of  each  contributor,  M^  and  N^  ,  can  readily  be  seen 
from  the  stress  curves  found  in  Figures  14  through  19.   Table  IV 
lists  the  maximum  non-dimensional  values  found  from  these  curves. 
The  bottom  of  the  ring,  ^  "  O  ,  is  seen  to  be  the  critical  stress 
location.   A  non-dimensional  value  of  -3.45  is  the  maximum  stress 
calculated  from  the  solution  for  the  original  ring.   In  comparison, 
the  modified  ring  solutions  have  values  of  -3.49,  -3.28,  and  -3.57, 
all  of  which  are  within  5  per  cent  of  the  original  ring  result. 
Consequently,  the  assumption  of  an  artificial  tangential  density 
has  a  negligible  effect  on  the  ring's  maximum  dynamic  response  for 
the  impulsive  loading  condition  considered. 
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In  conclusion,  the  proposed  method  for  reducing  the  stability 
restrictions  presented  in  Chapter  I  appears  to  be  valid.   The  same 
method  when  applied  to  the  finite  difference  approach  should  produce 
results  very  similar  to  those  presented  here.   Thus,  because  of  the 
assumed  artificial  tangential  inertia,  an  accurate  finite  difference 
solution  should  be  obtained  with  the  computation  time  reduced  by 
about  30  per  cent  for  ^  =  2,  40  per  cent  for  0  =  3,  and  50  per  cent 
for  ^=4  from  that  of  the  unmodified  ring  analysis. 
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APPENDIX  A 
SOLUTION  FOR  THE  SPECIFIC  IMPULSE  LOADING 
The  derivation  for  the  solution  to  the  radial  cosine  impulse 
is  outlined  in  this  appendix.   It  should  be  noted  that  the  develop- 
ment is  identical  to  that  of  the  exact  solution  when  the  factor 

y  =  1. 

Assuming  that  the  density  in  the  tangential  direction  is  larger 
by  a  factor  ^  than  that  associated  with  the  radial  direction,  the 
following  modification  in  the  ring  equations  will  result 


Ma<^^     -^j-  ^^^"^  ^T^ 


(A-l) 


Solutions  for  W  and  V  were  assumed  earlier  in  Equations  (11a) 
and  (lib)  and  are  repeated  below. 


(A-2) 


W(<^>r)  =  )     A„(t)  cos  n4> 


oo 


(A-3) 


v(4>  ,r)=  y   Q^{T)  sin  n4> 


The  symmetric  and  anti-symmetric  forms  for  W  and  V  arise  from  the 
symmetry  of  the  impulse  and  the  definition  of  the  V,  W  coordinate 
system. 
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Substituting  the  assumed  solutions  Into  the  modified  ring 
equations,  Equations  (A-1)  and  (A~2) ,  and  setting  the  two  algtbrale 
expressions  equal  to  zero  for  each  value  of  n  lead  to 

6T 


D(A,-r)Bj  +  A„-^(n^-O^AV^^  -o 


(A'5) 


for  n  "■  1 ,  2 ,  .  .  .  ,  oo 


Harmonic  forms  are  assumed  for  A  and  B 

n      n 


^„iT)^c.„J'^^ 


6„(T)  =  b„e 


(WbT 


(A- 6) 


where  a  and  b  are  arbitrary  constants.   Substitution  of  thti© 
n     n  "^ 

assumed  solutions  for  the  time  functions  into  Equations  (A-4)  and 

(A-5)  yields  the  following  two  expressions  from  which  the  common 

iodm  T 

time  term,  e   ^     ,  has  b©en  eliminated 


The  coefficients  a  and  b  cannot  be  determtned  from  the  above 

n     n 

equations  since  a  third  unknown,  Cd^  ,   is   involved,  Howevfr,  thty 
may  be  expressed  in  ternas  of  one  another  and  ellmlnattd  froBi  beth 
expressions. 
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From  Equation    (A-7) 


^n=   b, 


n  - 


Dn 


(A-9) 


Replacing  a   in  Equation  (A-8)  with  the  expression  above,  the 
n 

constant  b   is  common  to  all  terms  and  may  be  factored  out.   A 
n 

single  equation  in  the  unknown  CO   results 


n    CO., 


x,^°<(^^'0%n 


n 


n 


IcOr^^ 


+  nw  (n^-/)    =  0 


(A-IO) 


where 


^^=T 


r>? 


o<  = 


K 


bA^ 


The  modal  frequency  6l)   for  the  time  functions  A  and  B  may 

n  n      n 

now  be  determined.   Although  four  values  for  6i3   will  satisfy 
Equation  (A-10) ,  only  two,  the  positive  roots,  have  physical  signifi- 
cance.  These  roots  indicate  that  each  mode  of  vibration  is  character- 
ized by  two  distinct  frequencies 


COf)     =   -fy,  60c 


Cdn    =^  <5>n  CO, 


(A-ll) 


*  The  bars  indicate  that  the  frequency  coefficients  are 
functions  of  J  . 
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where 


(A-12) 


The  frequency  coefficients  F  and  g^^  may  be  simplified  In  light  of 
the  thin  ring  assumption.   For  small  values  of  C^  *,      these 
coefficients  are  approximated  by 


?n  -   ^(^'"Oi^ 


(A-13) 


The  initial  values  for  the  time  functions  must  now  be  determined 


in  order  to  solve  for  the  constants  &     and  b  .  At  time  T  ■  0.  thf 
expressions  for  the  displatements  and  veloeltlts  are  simple  Fourttr 
..series  expansions  In  the  single  variable  <^.   The  Fourier  coeffi- 
i clents  for  these  series  may  be  evaluated  to  yield  the  following 


initial  conditions  for  A^(T)  and  B  (T) 

n        n 

A  (0)  »  0 
n 


n  =»  0,  L ,  I ,    .  .  > 


*For  the  problem  considered,  o<*  3.54  x  10  . 
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<^A, 


Vo 


^(•^^--t 


'4^-1 


n  =  2,    4,    6, 


=  0 


n   =   3,    5,    7, 


B„(o)»0 


n  =   1 ,    2 ,    3 , 


#'»)■  ^ 


Mb  (0)  =  0 


n   =   2,    3,    4, 


(A-14) 
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The  time  functions  may  now  be  determined  from  the  above  initial 
values.   For  the  cases  n  =  0  and  n  =  1,  a  simple  substitution  into 
the  modified  ring  equations  yields 


Ao  (T)  =  — ^  S/A/  COoT 


(A-15) 


and 


^  ^  I  (A-16) 

The  solution  for  the  higher  mode  shapes,  n  >  1,  requires  further 
explanation.   Recalling  that  the  modal  frequency  (Oy^   was  character- 
ized by  four  distinct  roots,  there  are  four  separate  solutions  for 

A  and  B  which  satisfy  the  ring  equations.   The  sums  of  these 
n      n  on 

solutions  are  also  solutions.   Therefore,  the  general  time  functions 
for  n  >  1  are 

An (T)=  an,  e  "      -^0.^2  e  +  a^g  e  "*     *  -^  a,^  e 

(A-17) 


Although  eight  constants  appear  in  the  above  equations,  only  four 
are  independent.   As  determined  earlier  in  Equation  (A-9) ,  the  b 
constants  may  be  replaced  with  expressions  involving  the  "a  constants, 
Applying  Equations  (A- 13),  these  relationships  are 


58 


•^m  =  ^n, 


Y) 


n^-y^^ 


b^2-  ^H2 


r) 


!^"'  ^Cj 


\      (A-18) 


•^ns  "   ^"3 


n 


—  2 


n'-)rg^' 


b„4.=  a 


K)4 


n 


n^-ygrf 


Substituting  Equation  (A-18)  into  Equation  (A-17),  the  expressions 

for  A  (T)  and  B  (T)  contain  four  unknown  constants  to  which  four 
n         n 

initial  conditions  may  be  imposed. 

To  apply  the  initial  conditions,  Equations  (A-14),  the  time 
functions  must  be  expressed  in  terras  of  sines  and  cosines.   Making 
the  appropriate  expansions  for  the  exponentials  and  regrouping  the 
a  constants  into  new  constants,  a,  Equations  (A-17)  become 


Ar)  (T)  =    ^w;  C0$  fr^Cd^T  -H  ^r,z  S/V  ^n  ^^oT 


Bn(7-J=  -i^^   Qn,  COS  ?„0).T  ^  a^^  ^//V  f„  a),rj 


>(A-19) 


+ 


7frf^[^n3  C^S  Sn^cT  ^  ZT,^  SIN  <3«  ^ot] 


A  solution  may  now  be  determined  for  the  a  constants.   Applying  the 
initial  conditions,  the  cosine  terms  must  vanish.   Thus, 

5r„  =  3^3  =  b,,^   =   b„^  =  d?  (A-20) 

The  coefficients  for  the  sine  terms  are  given  by  the  following 
expressions 
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C^nz  = 


a 


104- 


2V.(-/f'"'(iff.^-K)^) 

)fw.7r(n^-0(f„^-S'f)-r„ 


n  =  2 ,  4,  6,  .  . 
(A-21) 


and  hence 


'r)Z 


'K)4. 


^r)2  =  ^n4.  =  ^ 


-2\/.(-/;^^'  n 


n  -  3,  5 ,  7 ,  . . 
(A-22) 


The  time  functions  may  now  be  written  directly  as 


n  =  2  ,  4 ,  6  .  .  . 
(A-23) 


n  =  3,  5,  7,  . . 
(A-24) 


An(T)  = 


2\4(:-/)^^' 


7:2 


(J-fi^-n' 


-2 


6„(T-)= 


even  modes  ;>  2 

At^  (T)  =  Br)  (t)  »  ^  °^^  "'^^^^  >  3 


(A-25) 
(A-26) 

These  expressions  may  be  substituted  into  Equation  (A-3)  to  obtain 
the  middle  surface  displacements,  W  and  V,  as  functions  of  time  and 
ring  position.   As  determined  earlier  in  Chapter  II,  Equations  (4) 
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and  (5) ,  the  stress-resultants  M  and  N  are  given  in  terms  of  the 
displacements  W  and  V  by 


N-^(i-;-w-i(ii.*w 


M  - 


K_l^'W 


i  (A-27) 


A^\^4>' 


■h  IV 


Thus,  the  solutions  for  W,  V,  M,  and  N  due  to  the  radial  cosine 
impulse  become 


W  %^    =     SIN  U>oT  - 


7T 


V 


^\[^ 


S/N\/Th~Jf(jO,T  COS   4, 


+ 


Z(-l) 


yi/2^1 


HYi'-m-^^) 


SIN  -f^  CdoT 


4-  = SIN    gn^oT 


COS    ncp  (A-28) 


V 


4fi7y^ 


SINp  -HJ^'  60«r    S/A/  cp 


oO 


+ 


zC~0 


%^l 


r)=z,4,-. 


Hn'-Di^^'-^-) 


-I- 


f) 


^  SIM  ^Yi  a),T 


S/A/   /^(^ 


(A-29) 
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and 


M  =-M 


k  IT  Ode 


oo 


=     SIN  cOoF  -^ 


n=2,-?,..- 


-2  (-/)"/^-' 


r    7:2 


l^n^  s;A/^;w„r 


-2 


3«  ^ 


f. 


^OS  n<?^ 


(A-30) 


Mo  =    A/^^f^^  =    -5/A/  U).T  -h    YW^   ^'^  \AH''^°^  ^^^  "^ 


A""      2  (-of 

L    (n^-/)(^7-Sv?) 


-■^  S//V  -^7^0 7" 


+  Sfvi   S//V  9kj  o^oT 


(T^S  ^<?b    -  o<  Mc 


(A-31) 
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APPENDIX   B 
CALCULATED   FREQUENCY  COEFFICIENTS 

TABLE   V 
FREQUENCY  COEFFICIENTS    (   ^  =   1) 


n 

f 
n 

Sn 

2 

0.016 

2.236 

4 

0.088 

4.123 

6 

0.208 

6.082 

8 

0.376 

8.062 

10 

0.592 

10.050 

12 

0.857 

12.041 

14 

1.169 

14.036 

16 

1.530 

16.031 

TABLE    VI 
FREQUENCY  COEFFICIENTS    (    g'  =   2) 


n 

f 
n 

^n 

2 

0. 

,015 

1. 

732 

4 

0. 

,085 

3, 

.000 

6 

0. 

,205 

4. 

,359 

8 

0. 

,373 

5. 

,745 

10 

0. 

,589 

7. 

,141 

12 

0. 

,854 

8. 

,544 

14 

1. 

,167 

9. 

,950 

16 

1. 

,528 

11. 

358 
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TABLE    VII 
FREQUENCY   COEFFICIENTS    (    U  =   3) 


n 

f 
n 

^n 

2 

0.014 

1.528 

4 

0.083 

2.517 

6 

0.202 

3.606 

8 

0.370 

4.726 

10 

0.587 

5.859 

12 

0.851 

7.000 

14 

1.163 

8.145 

16 

1.524 

9.291 

TABLE   VII 
FREQUENCY  COEFFICIENTS    (    g  =  4) 


n 

f 
n 

^n 

2 

0.013 

1.414 

4 

0.081 

2.236 

6 

0.199 

3.162 

8 

0.367 

4.123 

10 

0.583 

5.099 

12 

0.848 

6.082 

14 

1.161 

7.071 

16 

1.521 

8.062 
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